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ABSTRACT 

The three-dimensional (3D) structures of many 
biomacromolecules have been solved to reveal the 
functions of these molecules. However, these 3D 
structures have rarely been applied to constructing 
efficient molecular devices that function in living 
cells. Here, we demonstrate a 3D structure-based 
molecular design principle for constructing short 
hairpin RNA (shRNA)-mediated genetic information 
converters; these converters respond to specific 
proteins and trigger the desired gene expression 
by modulating the function of the RNA-processing 
enzyme Dicer. The inhibitory effect on Dicer 
cleavage against the shRNA designed to specifically 
bind to U1A spliceosomal protein was correlated 
with the degree of steric hindrance between Dicer 
and the shRNA-protein complex in vitro: The level of 
the hindrance was predicted based on the models. 
Moreover, the regulation of gene expression was 
achieved by using the shRNA converters designed 
to bind to the target U1A or nuclear factor-wB 
(NF-kB) p50 proteins expressed in human cells. 
The 3D molecular design approach is widely applic- 
able for developing new devices in synthetic 
biology. 

INTRODUCTION 

Artificial signal-converting systems have the potential to 
respond to intracellular biomacromolecules (e.g. specific 
proteins) and can control and rewire gene networks in 
living cells (l^t). These systems could be useful devices in 
biological science, medical science and engineering. Gener- 
ally, synthetic molecular devices have been designed by 
connecting naturally occurring modular parts [e.g. deoxy- 
ribonucleic acid (DNA), ribonucleic acid (RNA) and 



amino acid sequences] (5,6). Thus far, the method for 
constructing these devices depends on the primary and 
secondary structures of the biomolecules, without account- 
ing for their three-dimensional (3D) structures. For 
example, several synthetic RNA devices that control gene 
expression through the interaction between RNA and its 
ligands (e.g. antibiotics or proteins) have been designed 
and constructed (1,2,7-10). Specific target-sensing sites 
(e.g. naturally occurring RNA motifs or in vitro selected 
RNA aptamers) were inserted into RNA sequences based 
on the free energies of the RNA secondary structures 
(11-13). However, the capabilities of such devices are dif- 
ficult to predict in most cases. Another approach has been 
to use experimental evolution techniques to generate the 
functions of the RNA devices (14,15). These techniques 
are labor intensive in general, and it is impossible to 
predict whether the desired sequences are present in an 
RNA pool consisting of random sequences before evolu- 
tion. Consequently, it has been impossible to predict or 
optimize the performance of the devices without perform- 
ing a functional screen (14,15). 

As a new approach to constructing predictable and 
functional devices, 3D molecular design incorporating 
tertiary interactions between the biomacromolecular 
components of the device now seems to be feasible 
because of the increasing number of molecular com- 
plexes in the 3D structure-based library. 3D design of 
biomacromolecular complexes is likely to serve as an 
advanced method for producing next generation devices 
in synthetic and chemical biology because the perform- 
ance of such devices is highly predictable. However, a 
rational 3D design has not yet been implemented for con- 
structing useful biomacromolecular complexes [e.g. RNA 
protein complexes (RNP)] that function in living cells 
(11,16,17). 

Naturally occurring RNA-mediated signal converters 
have been discovered in human systems (18-21). For 
example, the Lin28 protein interacts with pre-/e/-7 
microRNA (miRNA) and blocks the processing of 
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pre-let-7 miRNA into mature let-7 miRNA, which leads to 
the derepression of let-7 target genes, including oncogenes 
and cell cycle genes. Crystal and nuclear magnetic reson- 
ance structural analyses have shown that the cooperative 
interactions of two distinct regions in the Y\n2i-let-7 
complex are critical for inhibiting the activity of the 
RNA-processing enzyme Dicer (22). Accordingly, by 
mimicking the function of natural RNAs, we sought to con- 
struct a synthetic short hairpin RNA (shRNA) that 
responds to a protein of interest and activates target gene 
expression. 

Here, we report the 3D structure-based molecular design 
and construction of RNA-mediated protein information 
converters. The RNA device is designed to interact with 
a target protein that is expressed in cells and trigger the 
desired gene expression. We designed protein-responsive 
shRNAs in silico on the basis of high-resolution structures 
of known RNP complexes and evaluated their perform- 
ance in mammalian cells. The designed shRNA senses 
the target proteins and inhibits the activity of Dicer, 
thereby controlling the RNA interference (RNAi) 
pathway in a predictable manner. 

MATERIALS AND METHODS 

3D modeling of protein responsive shRNAs and their 
trigger proteins 

The 3D atomic model of the N-terminal domain of U1A 
(human) and its binding motif in UlsnRNA (UlAi motif) 
or in the 3'-untranslated region of U1A messenger 
RNA (mRNA; U1A 2 motif) were obtained from Protein 
Data Bank (PDB) (ID: 1URN and ID: 1AUD, respect- 
ively). The 3D atomic model of p50 [mouse nuclear 
factor-KB (NF-kB)] bound to its aptamer RNA was 
obtained from PDB (ID: lOOA). The 3D model of 
human NF-kB p50 and its aptamer has yet to be 
determined; therefore, we used the model of mouse 
NF-kB p50 for our 3D design. Beforehand, we verified 
the amino acid sequence differences between mouse and 
human p50 and confirmed that the six point mutations 
between the two were on either the terminal region or the 
loop far from the RNA-binding domain; therefore, the mu- 
tations are likely not critical for the structure and binding 
affinity. Molecular designs were performed with Discovery 
Studio (Accelrys) as previously described (23). Parts of the 
protein sequences were trimmed in silico to correspond to 
the regions of p50 and Ul A-expressing proteins. The RNA 
stem region of the p50 aptamer and UlA-binding motifs 
were connected to double-stranded RNA (dsRNA) by 
using the superimpose protocol. The corresponding 
atomic coordinates of sugar, base and phosphate backbone 
in the four terminal bp of dsRNA and the protein-binding 
RNA motif were tethered, and the distance between these 
tethered coordinates was minimized using least squares ap- 
proximation polynomial. 

shRNA transcription in vitro 

For in vitro shRNA transcription, a single-stranded tem- 
plate was annealed to the T7 annealing primer 5'-GCTAA 



TACG ACTC ACTATA-3' . shRNAs were transcribed and 
purified as follows: the transcription reaction mixtures 
were incubated for 3-4 h at 37°C. Then, the template 
DNA was degraded for 20min at 37°C with 10 ul of 
TURBO DNase (Ambion). Any proteins in the mixture 
were removed by phenol-chloroform extraction, and nu- 
cleotide monomers were removed by PD-10 column puri- 
fication (GE Healthcare) and ethanol precipitation. The 
shRNAs were purified using 15% native polyacrylamide 
gel electrophoresis and eluted with 500 ul of elution buffer 
(0.5 M NaCl, 0.1% sodium dodecyl sulphate and ImM 
EDTA) at 37°C overnight. The eluted shRNAs were 
filtered with a 22- um microfilter (Millex GP) and 
precipitated with ethanol. The DNA templates used in 
this study are described in Supplementary Table S2. 

Inhibition assay of human Dicer cleavage 
for shRNA in vitro 

In vz'fro-transcribed U1A protein-responsive shRNA 
device variants (UlA r or dUlA r sh 21-31; 0.5 uM) and 
recombinant U1A protein (0 or 5uM) were mixed with 
1 ul of lOmM adenosine triphosphate (ATP), 0.5 ul of 
50 mM MgCl 2 , 4ul of Dicer Reaction Buffer (Genlantis), 
1 ul of 0.5 unit/ul Recombinant Human Dicer Enzyme 
(Genlantis) and water (up to 10 ul). The mixtures were 
incubated at 37°C for 15 h, and the incubation was 
stopped using 1 ul of Dicer stop solution. The mixtures 
were separated by 15% native polyacrylamide gel electro- 
phoresis at 4°C for 30min, and the gel was stained with 
SYBER-Green I and II (Takara Bio Inc). The fluorescence 
of the stained RNA was detected by FLA-7000 (Fujifilm). 
The assays using other samples were conducted using the 
above mentioned procedure. 

Construction of NF-kB (p50), U1A or MS2 coat protein 
expression plasmids 

The NF-kB p50 (amino acids 39-352 of human NF-kB1) 
gene was amplified by polymerase chain reaction from 
total complementary DNA (cDNA) from HeLa cells, 
using the following primers: NF-kB (p50) Fwd (5'-CCA 
GGATCCCACCATGGATGGCCCATACCTTCAAAT 
ATTAGAGC-3') and NF-kB (p50) Rev (5'-tcgCTCGAG 
TTCAGGATAGTAGAGGAAAGGTTTTGG-3). The 
gene was then cloned into pcDNA3.1-mycHisA 
(pcDNA, Invitrogen) between the BamHI and Xhol 
sites to construct the p50 expression plasmid, and the 
sequence was verified. 

The sequence coding U1A (amino acids 2-98 of human 
U1A, without nuclear localization signal) or MS2 (bac- 
teriophage MS2 coat protein with amino acids 68-81 
deleted) was amplified using the following primers: 
U1A_T Fwd (5'-GACACCATGGGATCCGCAGTTCC 
CGAGA -3')/UlA_T Rev (5'-TATGAGATCTTTTCAT 
CTTGGCAATGATATCTG-3') or MS2-CP Fwd (5'-CA 
CCATGGGATCCGCTTCTAACTTTACTCAGTTCGT 
TCTC-3')/MS2-CP Rev ( 5' -TATG AG ATCTGT AGATG 
CCGGAGTTGGC-3') from human U1A or MS2-GFP 
expression plasmids, respectively. In short, the U1A or 
MS2 expression plasmid was cloned into pcDNA5- 
mycHisA (pcDNA, Invitrogen) containing an Internal 



Nucleic Acids Research, 2012, Vol. 40, No. 18 9371 



Ribosomal Entry Site-(IRES)-driven DsRed-Express ex- 
pression cassette from pIRES2-DsRed-Express (Clontech 
Laboratories, Mountain View, CA, USA) to construct the 
U1A or MS2 expression plasmid, and the sequence was 
verified. 

EGFP expression control of the devices in 
cultured cells 

The oligo sequences for constructing the shRNA expres- 
sion plasmids are shown in (Supplementary Table SI). To 
assess the p50-responsive shRNAs (p50- or dp50-sh 22-28), 
293FT cells (1.5 x 10 5 per well) grown for one day in 
24-well plates were co-transfected with 0.3 ug of one of 
the following: p50- or dp50-sh 22-28 expression plasmid, 
0.3 ug p50 or MS2 expression plasmid and Enhanced Green 
Fluorescent Protein (EGFP) expression plasmid using 2 ul 
of Lipofectamine 2000 (Invitrogen). To assess the 
UlA-responsive shRNAs, we used the corresponding 
shRNA-expressing plasmid instead of p50-sh 22-28. 
Twenty-four hours after transfection, following fluorescent 
microscopic analysis, the cells were collected and resus- 
pended in the medium, and the mean intensity of EGFP 
fluorescence was determined using a flow cytometer 
(FACSAria, BD Bioscience). To exclude abnormal cells, 
those cells whose Foward Scatter and Side Scatter signal 
scopes were identical to those of untreated 293FT cells were 
selected for the calculation of mean intensities, and 30 000 
cells per each sample were analysed. To measure the EGFP 
signal, a 488-nm semiconductor laser beam and a 530/30- 
nm emission detector were used. 

RESULTS 

Design principle of protein-responsive shRNA devices 

We chose to design shRNA structures because they 
contain a simple stem-loop structure and have been used 
as powerful tools to knock down desired genes through 
the RNAi pathway. Human Dicer recognizes shRNA at 
the end of its stem and cleaves it at approximately the 
22nd nucleotide from the 3' and/or 5' end to produce a 
functional dsRNA for performing RNAi (24,25). Human 
Dicer does not discriminate according to the size and 
structure of the shRNA loop region during its processing; 
therefore, desired sequences can generally be incorporated 
into the loop region (26). We have previously constructed 
a synthetic shRNA by simply inserting a kink-turn RNA 
motif (Kt), which binds to the ribosomal protein L7Ae in 
Archaeoglobus fulgidus (2), into the shRNA loop region. 
L7Ae modulates the function of Dicer through its specific 
binding to the shRNA; however, it is unclear whether the 
protein-binding site is optimally set in the shRNA to 
maximize the steric hindrance to obstruct Dicer. To 
solve this problem, we took a 3D molecular design 
approach to constructing sophisticated protein-responsive 
shRNA devices. The ability of the RNA-binding protein 
to inhibit Dicer-mediated cleavage should directly deter- 
mine the RNA device's efficiency of RNAi control because 
steric hindrance can be predicted from the 3D model of 
Dicer and the shRNA-protein complexes (schematic rep- 
resentation and corresponding 3D structure are shown in 



Figure 1A). Accordingly, the engineered shRNA with a 
specific protein-binding site was designed by adding bp 
between the double-stranded region and the loop to 
adjust the orientation and distance from the Dicer 
cleavage site (Figure IB). 

In the 3D model, we predicted the region that is 
required by human Dicer to cleave the shRNA based on 
the high-resolution X-ray structure of Giardia Dicer and a 
single-particle analysis of human Dicer using electron mi- 
croscopy (27,28). Although the entire crystal structure of 
human Dicer has not been solved, previous structural 
studies have revealed that the regions responsible for 
nuclease activity of both human and Giardia Dicers are 
similar (27,29). Two structural models of human Dicer 
have been proposed by incorporating the Giardia Dicer 
crystal structure into the human Dicer electron micros- 
copy map (27). A credible model from previous structural 
and biochemical studies indicates that a helicase domain is 
located adjacent to the RNaselll domains (26,27). 
However, it has been reported that the disruption of the 
helicase domain maintains the stability and activity of 
human Dicer (30-32), indicating that the activity of 
human Dicer is not affected. Therefore, we did not 
account for the helicase domain in our 3D model. 

3D molecular design of UlA-responsive shRNA devices 

We designed a series of shRNA devices that respond to the 
expression of U1A (referred to as UlA r sh 21-31; Figure 
2). U1A is one of the most abundantly expressed 
mRNA-processing and housekeeping proteins in mamma- 
lian cells (33). We used the U1A N-terminal domain (amino 
acids 2-98), which does not contain the region for nuclear 
localization signal (amino acids 100-144), indicating that 
the expressed U1A will be present in both the cytoplasm 
and nucleus. The UlA r sh 21-31 RNAs consist of a loop 
region containing an RNA motif for binding to U1A and 
21-31 consecutive bp containing the target sequence for an 
mRNA-encoding EGFP (Supplementary Figure S1A). The 
engineered shRNAs were designed based on a specific RNP 
complex, which is UlA's N-terminal domain that is 
associated with the corresponding UlA-binding motif in 
human Ul's small nuclear RNA stem-loop II, whose 
high-resolution 3D structure is known (34,35). We used 
a 3D modeling method to define the degree of steric 
hindrance between human Dicer and U1A to optimize 
the performance of the device (+U1A, Figure 2, 
Supplementary Movies SI and S2). Thus, if U1A interfered 
with the region that is required for the Dicer's recognition, 
it is likely to cause inhibition of the Dicer function. Indeed, 
the steric hindrance between the Dicer and U1A (with 
UlA r sh 23-29) was observed, and then, the RNP 
complex of UlAi-sh 30-31 exceeded the expected access 
region of Dicer in our in silico 3D modeling (Figure 2, 
Supplementary Movie S2). 

Biochemical and in cell analyses of UlA-responsive 
shRNA devices 

To confirm the prediction that was derived from the 3D 
modeling of the UlA r shRNA device, we investigated 
whether the in v/fro-transcribed UlA r sh 21-31 and the 
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Figure 1. Schematic and structural model of the protein-responsive shRNA device mechanism and its bound protein orientation. (A) In the absence 
of the RNA-binding protein, the shRNA-processing enzyme Dicer (grey) can access and process the designed shRNA, and the processed shRNA 
induces the knockdown of its target gene through an RNAi mechanism (left). In contrast, Dicer is inaccessible to shRNA in the presence of a protein 
(pink) that binds to the loop region of the shRNA (yellow) because the RNP interaction faces Dicer and inhibits its access. The prevention of Dicer's 
function causes the derepression of the gene knockdown (right). (B) Protein-responsive shRNA devices A and B are designed by connecting an 
RNA-protein interaction motif with each length (A: 24 bp, B: 21 bp) of dsRNA. By specifying the position of two devices by reference to Dicer 
cleavage sites and comparing these two bound proteins' locations, we found that the protein bound on device A is located ~7.8A more distant to the 
Dicer cleavage sites and rotates ~90° in an anti-clockwise direction around the axis of the dsRNA compared with that bound on device B. Dicer 
cleavage sites are shown in two filled red and blue circles. 



corresponding negative control, dUlA r sh 21-31, which 
contained a defective U1A binding motif (Supplementary 
Figure SIB), are processed by recombinant human Dicer 
in the presence or absence of U1A protein (Figure 3 A and 
B and Supplementary Figure S2). In the absence of U1A, 



Dicer cleaved the transcribed UlA r sh 21-31 in vitro, 
which indicated that human Dicer can recognize and 
cleave the designed shRNAs containing the UlA-binding 
motif at the loop region (Figure 3A). As expected, the 
predicted effect of steric hindrance corresponded to the 
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Figure 2. Predicted steric hindrance of U1A with Dicer in the 3D model of UlA-responsive shRNA devices. A 3D model of the catalytic domain 
and peripheral structure of Giardia Dicer (PDB ID: 2QVW, 326-385, 401- C-terminal amino acids), which is analogous to that of human Dicer, 
(white and white translucent) is displayed. Eleven U1A protein-responsive shRNA device variants, UlA r sh 21-31 (white ribbon) and bound U1A 
(orange) are shown. UlA r sh 21-31 consists of a loop region (green) that contains a UlA-binding motif with 21-31 bp of dsRNA targeting the 
EGFP gene. These models were arrayed with reference to the dsRNA axis and by producing two nucleotides (red spheres: 22nd nucleotide from the 
5' end and blue spheres: 22nd nucleotide from the 3' end) that are located between the two metal ions of catalytic sites of Giardia Dicer. The expected 
accessing region of Dicer was indicated by a white translucent area. 



degree of inhibition of the Dicer function in the presence 
of U1A (Figure 3 A and B and Supplementary Figure 
S2B). Importantly, the UlA r sh 30-31 less efficiently in- 
hibited Dicer cleavage when compared with UlA r sh 
24-29, indicating that the longer stem in UlA r sh 30-31 
can release the UlA-mediated inhibition of Dicer activity 
(Figure 3A and B). As anticipated, the result is consistent 
with the prediction based on the in silico models (Figure 2, 
sh30-31). We next investigated the cellular function of 
UlAj-sh 21-25 because the inhibition efficiency of Dicer 
cleavage observed for UlA r sh 26-29 stays roughly 
constant. The designed UlA r sh 21-25 expression 
plasmids or the corresponding negative control dUlAj- 
sh 21-25 expression plasmids were transfected into 
human cultured 293FT cells expressing U1A or the 
control bacteriophage MS2 coat protein. Twenty-four 
hours after transfection, the level of EGFP was 



measured with a flow cytometer. In the MS2-expressing 
cells, efficient down-regulation of EGFP expression was 
observed when UlA r sh 22-25 was transfected (red bars, 
Figure 3C). RNAi activity of UlA r sh 21 and dUlA r sh 
21 was weaker than that of the other constructs (sh 22-25) 
(Figure 3C), indicating that the 21 -bp length of the 
double-stranded region of the shRNAs is likely incompe- 
tent for Dicer processing in cell. The effective 
de-repression of EGFP expression when UlA r sh 25 
was transfected into UlA-expressing cells was observed 
by flow cytometry (blue bars, Figure 3C and 
Supplementary Figure S3) and fluorescent microscopy 
(Figure 3D). The binding-defective shRNA expression 
plasmids (dUlA r sh 22-25) had no influence on RNAi 
activity in vitro Dicer cleavage assay (Supplementary 
Figure S2) and in cells expressing U1A (Figure 3C 
and Supplementary Figure S3), indicating that the 



9374 Nucleic Acids Research, 2012, Vol. 40, No. 18 



A Dicer + + + + + + + + + + 

U1A(5mM) - + - + - + - + -+ -+ +. + . + .+ 



+ . + 



RNP . 

complex 



=3 



Cleaved 
stem loop 
+U1A 



shRNA 

Cleaved 
dsRNA . 



.unyyyy yuyutL 

— ^ w ^ M w I 

n * 



I 



Cleaved 
stem loop 



0.5pM ' 
U1A-sh 



B 




29 



30 



31 



31 



normalized band intensity of cleaved dsRNA 
[(U1A+) U1A,-sh / (U1A-) U1A r sh] / 
[(U1A+) dU1A,-sh / (U1A-) dU1A,-sh] 




0.5|jM U1A-sh21 



22 23 24 25 26 27 28 29 30 31 



■ U1A + U1A,-sh .MS2 + U1A,-sh 
"U1A + dU1A,-sh "MS2 + dU1A,-sh 





U1A + dU1A,-sh25 



MS2 + dU1A-sh25 



Figure 3. In vitro and cellular performance of UlA-responsive shRNA devices. (A) Human Dicer cleavage assay in vitro. UlA r sh 21-31 and U1A 
formed the RNP complex. UlA^sh 23-29 modulated Dicer activity in the presence of U1A. (B) The relative efficiency values of Dicer cleavage 
observed in (A) were calculated using the ratio of the cleaved band intensity in the presence of U1A to that in the absence of U1A and normalized to 
the value of dUlA^sh 21-31 (see also Supplementary Figure S2). (C) Flow cytometric analysis of the regulation of the RNA device variants by U1A 
in 293FT cells. 293FT cells were co-transfected with the three plasmids that expressed U1A or MS2, EGFP and one of the U1A|- or dUlA r sh 
variants and were analysed the next day. The relative expression levels of EGFP fluorescence in the U1A- or MS2-expressing cells were normalized to 
the EGFP expression levels of the sh-N (negative control for EGFP knockdown, see Supplementary Table SI) and U1A- or MS2-expressing cells. 
The error bars indicate the standard deviations of three independent experiments. sh-GFP was used as a positive control for EGFP knockdown. (D) 
Fluorescent microscopic images of cells that expressed EGFP, UlA r or dU!A r sh 25 and U1A or MS2 from each plasmid. 



RNA-protein interaction is responsible for suppressing 
RNAi activity. Thus, it is conceivable that the degree of 
steric hindrance between U1A and human Dicer is pre- 
dictable from the 3D models. U1A located on the UlA r 
sh 25 effectively interfered with the Dicer accessible region 
(Figure 2 and Supplementary Movie S2). The effect of this 
steric hindrance observed in UlA r sh23-25-UlA complex 
models was positively correlated with the ability to 



de-repress EGFP expression in a phased manner 
(Figures 2 and 3 and Supplementary Figure S3). 

3D molecular design and in cell analyses of p50-responsive 
shRNA devices 

To confirm the versatility of the molecular design based on 
the 3D structure, we constructed a series of shRNA 
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devices that respond to the p50 domain of NF-kB (Figure 
4A and Supplementary Movies S3 and S4). NF-kB is one 
of the major transcription factors in eukaryotic cells, and 
its abnormal expression has been linked to the prolifer- 
ation of several types of tumor cells (36,37). The 
NF-KB-responsive shRNAs (referred to as p50-sh 22-28) 
consist of a loop region containing a specific binding site 
for p50 and a stem (22-28 bp), which contains the target- 
ing sequence for an mRNA-coding EGFP (Figure 4A and 
Supplementary Figure S4A). The p50 domain consists of 
two parts: the N-terminal domain (p50-N, colored sky 
blue) and the C-terminal domain (p50-C, colored green; 
p50-N and p50-C are connected via a linker peptide in 
Figure 4A). The p50-N, but not the p50-C, domain 
tightly binds to the aptamer sequence inserted into the 
loop region of the shRNA (Figure 4A) (38,39). Notably, 
p50 does not contain a nuclear localization signal peptide; 
therefore, the expressed p50 will be present in the both the 
cytoplasm and nucleus. 

The properties of the designed p50-sh 22-28 expres- 
sion plasmids and the corresponding negative control 
dp50-sh 22-28 expression plasmids containing a defective 
p50 aptamer (Supplementary Figure S4B) were examined 
in 293FT cells that expressed p50 or MS2. In the 
MS2-expressing cells, a pronounced and moderate 

A p50-sh 22 23 24 



down-regulation of EGFP expression was observed 
when p50-sh 22-26 and p50-sh 27-28 expression 
plasmids were transfected, respectively (Figure 4B). 
In vz/ro-transcribed p50-sh 22-28 RNAs are processed 
by recombinant human Dicer in a similar manner 
(Supplementary Figure S5A and S5B). In the presence of 
p50, the down-regulation of EGFP was suppressed 
when p50-sh 26-28 expression plasmids were transfected 
(Figure 4B and C and Supplementary Figure S5C). The 
defective shRNAs (dp50-sh 22-28) that were incapable of 
binding to p50 had no influence on RNAi activity in cells 
expressing p50 (Figure 4B & Supplementary Figure S5C). 
It is conceivable that the degree of the steric hindrance 
between p50 and Dicer is predictable from the 3D 
models as in the case of U1A; p50 that is located on sh 
26-28 effectively interferes with the Dicer accessible region 
(+p50, Figure 4A and Supplementary Movie S4). 

3D design and cellular function of an alternative device 
with the protein-stabilized RNA structural motif 

As an extension of this design, we constructed an alterna- 
tive device with a more complex 3D RNP structure 
(Supplementary Figure S6). We attempted to inhibit 
Dicer function by using both an RNA-binding protein 
and the protein-stabilized RNA structural motif because 

25 26 27 28 
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p50(dp50)-sh 



p50 + dp50-sh27 MS2 + dp50-sh27 



Figure 4. Prediction and assessment of the steric hindrance between Dicer and NF-kB p50 attached to p50-responsive shRNA devices. (A) The 3D 
model of seven p50 protein-responsive shRNA device variants, p50-sh 22-28 (white ribbon) and bound p50 (N-terminal domain with linker: sky blue; 
C-terminal domain: dark green). p50-sh 22-28 consists of the loop region (yellow) that contains an aptamer for p50 with 22-28 bp of dsRNA. These 
models were arrayed using the same procedure as shown in Figure 2. (B) Flow cytometric analysis of the regulation of the p50-sh variants by p50 
expression in 293FT cells. 293FT cells were co-transfected with the three plasmids that expressed p50 or MS2, EGFP and one of the p50- or dp50-sh 
variants and were analysed the next day. The relative expression levels of EGFP fluorescence in the p50- or MS2-expressing cells were normalized to 
the EGFP expression levels of the sh-N and p50 or MS2 expressing cells. The error bars indicate the standard deviations of three independent 
experiments. sh-GFP was used as the positive control for EGFP knockdown. (C) Fluorescent microscopic images of cells that expressed EGFP, 
p50- or dp50-sh 27, and p50 or MS2 from each plasmid. 
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a certain RNA-protein interaction is known to stabilize 
the RNA structural motif containing the protein-binding 
site (40). Accordingly, we incorporated the UlA-binding 
RNA motif into the loop region that is flanked by the two 
RNA stems; the UlA-stabilized RNA structural motif 
presumably contributes to the effect of steric hindrance. 
The corresponding construct (UlA 2 -sh 21-27) was 
designed to exhibit a loop-stem-loop structure (shown in 
green) and contains the UlA-binding motif that is con- 
nected to the main stem (21-27 bp; Supplementary Figure 
S6A) (33,40). The loop-stem-loop motif of UlA 2 -sh is pre- 
dicted to be more rigid than the loop structure of Ul A r sh 
in terms of RNA secondary structure (40). The 3D model 
of the shRNA-UlA complex shows that the protruding 
portion of the RNA UUCG loop with the 5'-GUCC/ 
GGAC-3' stem region (which is referred to as the 
UUCG stem-loop) and U1A are arranged symmetrically 
(Supplementary Figure S6B). The steric hindrance 
between Dicer and U1A (UlA 2 -sh 21 and UlA 2 -sh 
24-27) or between Dicer and the UUCG stem-loop that 
is bound to U1A (UlA 2 -sh 21-24) can be seen in the 3D 
models. 

To test the cellular function of the UlA 2 -sh 21-27 RNA 
devices, we transfected the designed expression plasmids 
(UlA 2 -sh 21-27 and dUlA 2 -sh21-27; Supplementary 
Figure S6A) into 293FT cells that expressed either U1A 
or control MS2 (Supplementary Figure S7). The 
suppression of RNAi was observed when the UlA 2 -sh 
22-26 constructs were transfected into the 
UlA-expressing cells. We confirmed that the designed 
UlA 2 -sh 24 interacted specifically with recombinant 
U1A and that the catalytic activity of Dicer was efficiently 
blocked by the RNP interaction in vitro (Supplementary 
Figure S8). U1A presumably stabilizes the conformation 
of the UUCG stem-loop region so that both the 
UlA-stabilized UUCG stem-loop region and U1A 
should become responsible for steric hindrance with 
Dicer (Supplementary Figure S6). Both UlA 2 -sh 21 and 
UlA 2 -sh 27 and their negative control plasmids induce 
RNAi weakly in the absence of U1A when compared 
with UlA 2 -sh 22-26 and their negative control plasmids, 
indicating that the performance of the designed switches 
reflect a balance between the proper positioning of the 
UlA-shRNA complex to interfere with Dicer and the 
activity of the designed shRNA (without protein) to 
proceed in the correct RNAi pathway. In particular, 
both UlA 2 -sh 27 and p50-sh 27-28 (and their negative 
control plasmids) with longer dsRNA region induce 
RNAi weakly in the absence of the corresponding 
protein compared with the other constructs with shorter 
dsRNA regions (Figure 4B and Supplementary Figure 
S7A). Thus, longer dsRNA region (>27bp) of the 
shRNA seems incompetent for facilitating cellular RNAi 
pathway. Consequently, synthetic shRNAs with the 
proper lengths of dsRNA (22-27 bp) should be used for 
building a more diverse array of orthogonal devices: the 
range from 22 to 27 bp dsRNA is likely sufficient to 
generate the optimal shRNA devices. The 22-27 
shRNAs can rotate the RNA motif -150° (Figures 2 
and 4A) so that the steric hindrance can be optimized 
within the range in most cases. 



DISCUSSION 

We have demonstrated that the 3D modeling of synthetic 
shRNAs and RNA-binding proteins is highly useful for 
developing efficient and tunable protein-responsive RNA 
devices; in silico 3D modeling predicted functional devices 
that control RNAi in mammalian cells. The designed 
protein-responsive shRNAs provide the proper orienta- 
tion of the bound protein on the structured RNA loop 
for disturbing Dicer-dependent RNA cleavage activity in 
the 3D model. 

For our 3D model of human Dicer, we referred to 
RNaselU's catalytic domain and the peripheral structure 
of Giardia Dicer because structural studies have indicated 
that the crystal structure of Giardia Dicer is geometrically 
similar to the nuclease core of human Dicer (27,29). For 
our in silico 3D modeling, in vitro Dicer cleavage and 
cellular gene expression analyses indicated that there is a 
clear correlation between the model and our experimental 
data: The designed RNA devices can efficiently regulate 
the function of Dicer in human cells as in the case of the 
in vitro experiments, although we are unable to disregard 
the possibility that unexpected interactions between the 
RNP and some other molecules might affect the regula- 
tion of the shRNA devices in cells. 

The introduction of a set of RNA devices that respond 
to multiple proteins could simultaneously and inde- 
pendently control multiple signals to improve the per- 
formance of synthetic signaling circuits in a single cell. 
For example, we have succeeded in re-directing apoptosis 
pathways of cultured human cells by using two distinct 
RNA devices that simultaneously regulate the expression 
of two apoptosis regulatory proteins (2). In addition, the 
RNA devices responding to endogenously produced 
proteins could rewire gene regulatory networks for thera- 
peutic applications. For example, the in vivo systemic 
delivery systems of shRNA have been developed, 
although the precise delivery to the specific cells of 
interest remains unresolved (41). Our protein-responsive 
shRNA devices could be advantageous for controlling 
RNAi activity in the target cells because the function of 
synthetic shRNAs can be regulated depending on intracel- 
lular environment. 

Precisely and strictly controllable RNA devices could be 
designed and developed based on the 3D structure by: 
(i) introducing multiple RNP interaction motifs into the 
appropriate region(s) of the shRNA device, (ii) altering 
the binding affinity and specificity to the sensing protein 
of interest or (hi) regulating the degree of steric hindrance 
between RNA and the corresponding RNA-binding 
proteins. Therefore, the 3D design of synthetic RNP- 
mediated information converters could contribute to the 
synthetic biology and other fields as a versatile tool for 
biomacromolecular design (9,16,17). This design principle 
can be applied to other biomacromolecular complexes 
(e.g. DNA-protein or protein-protein complexes) whose 
high-resolution 3D structures have already been solved. In 
addition, the 3D design strategies that are used to inves- 
tigate the effects of steric hindrance between RNA and 
proteins on the function of RNA-processing enzymes 
could be useful for studying the relationship between the 
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molecular mechanisms and the structures of naturally 
occurring RNP complexes that are involved in RNA pro- 
cessing and/or regulation pathways (22,27). 
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Supplementary Data are available at NAR Online: 
Supplementary Tables 1 and 2, Supplementary Figures 
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